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Abstract. We study a fermionic version of the Sherrington-Kirkpatrick model including nearest-neighbor 
hopping on a oo-dimensional simple cubic lattices. The problem is reduced to one of free fermions moving 
in a dynamical effective random medium. By means of a CPA method we derive a set of self-consistency 
equations for the spin glass order parameter and for the Fourier components of the local spin susceptibility. 
In order to solve these equations numerically we employ an approximation scheme which restricts the 
dynamics to a feasible number of the leading Fourier components. From a sequence of systematically 
improved dynamical approximations we estimate the location of the quantum critical point. 
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1 Introduction 

Quantum-dynamical mean field theories are highly ap- 
preciated, helpful methods to provide insight into models 
of strongly interacting systems 1 . A well-known case is 
certainly the Hubbard model including its various model 
ramifications. Fermionic spin glass models, though describ- 
ing the quite different physics of randomly interacting dis- 
ordered systems, belong to this category as well; in par- 
ticular, they also deal with an intense interplay between 
magnetism and transport properties. 

Quantum dynamics in systems with random many- 
body interactions are known to contain a challenging tech- 
nical difficulty which arises from the double-time depen- 
dence of the spin glass field. This feature of quantum spin 
dynamics is hard to handle even in the dynamical mean 
field theory (DMFT) and requires further simplifications, 
such as the limit M — > oo in S'C/(M)-generalizations of 
spin systems 

In earlier work the so-called static approximation (de- 
fined by a static effective spin glass field) was used to con- 
struct a systematic low temperature expansion for metal- 
lic spin glasses j^]- The quantum critical point was de- 
termined within this approximation. Later, a Ginzburg- 
Landau theory 0] , which kept only terms relevant or dan- 
gerously irrelevant under the quantum dynamical renor- 
malization group, was applied to study the critical expo- 
nents. 

Within this framework it turned out that the quantum- 
dynamical Gaussian result for the shift exponent of the 
critical temperature disagreed with the static approxima- 
tion result due to the effect of dangerously irrelevant dy- 
namic couplings. The shift exponent 4> can be defined by 



T c (x) ~ (x c - x)i 



> x, where x stands for the pa- 



rameter which drives the quantum phase transition, and 
T c {x c ) vanishes by definition at the quantum critical point 
(QCP). 

Naturally, renormalization groups are not designed to 
determine critical points. However, the positions of the 
QCPs are relevant, too. In finite dimensions the magnetic 
transition may coincide with a localization transition driv- 
en by the randomness of the many-body interaction. If 
these transitions do not coincide it is still important to 
know whether or not the magnetic transition occurs within 
the metallic phase and, perhaps, to identify parameters 
which influence the relative position of the QCPs. The 
answer may in some detail be model-dependent, but a co- 
incidence of both transitions would almost certainly con- 
tain deeper reasons and could be expected to appear at 
least in model classes. 

In this work we adapt and apply the coherent poten- 
tial approximation (CPA) to a metallic spin glass problem. 
Originally this powerful non-perturbative method was de- 
veloped to describe non-interacting disordered electron 
systems [5], and in this context the CPA can be shown 
to become exact in the limit of infinite spatial dimen- 
sions (d — > oo) Later, the CPA formalism was gener- 
alized to deal with interacting electron systems |7I8| with 
highly non-trivial couplings of the Matsubara frequencies. 
Finally, the CPA method can also be applied to dynam- 
ical disorder. Like in the present work, this situation for 
instance results from the dynamical decoupling of inter- 
action terms 

It is interesting to note that effective action terms 
that are non-diagonal in imaginary time not only occur 
in the presence of disorder. For certain translationally in- 
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variant models, the so-called extended DMFT approach 
[TO] , which to some extend incorporates spatial correla- 
tions, gives rise to effective impurity problems equivalent 
to the one studied in this article. 

The present article is organized as follows. The re- 
quired technical tools and strategies are developed in Sec. 
After the model definition we give a brief description 
of the dynamical two-step decoupling procedure, which 
reduces the problem to one of free fermions moving in 
a dynamical random magnetic field. The model assump- 
tion of a fully connected magnetic interaction facilitates 
a saddle point treatment. Extending previous work |11I3| 
we choose a dynamical self-consistent saddle point for the 
spin glass field. 

The quantum-dynamical CPA method, shaped to ap- 
ply to the present quantum spin glass model, is discussed 
in Sec. 12. 21 Within this framework, a set of replica-symmet- 
ric self-consistency equations for the saddle point values 
q (spin glass order parameter) and q m = q(ui m ) (the local 
replica diagonal spin correlation) are derived. As the cen- 
terpiece we obtain a matrix CPA-equation in Matsubara 
frequency space. 

Our ideas for the approximate solution of these self- 
consistency equations are introduced in Sec. 12.31 As a sys- 
tematic approximation scheme we propose to restrict the 
dynamics of the effective random medium to a number 
of bosonic Matsubara frequencies that can be numerically 
dealt with. Similar approximations have been constructed 
earlier in the context of the Ising spin glass in the trans- 
verse field by means of discretization of the imaginary time 
axis [HEj. 

Section |3 revisits the spin-static approximation where 
a crucial simplification arises: the matrix structure of the 
equations disappears. In particular, the CPA-equation can 
be solved independently for each Matsubara frequency. 
However, via the saddle point values, the non-trivial cou- 
pling of the frequencies is preserved also in the static set 
of self-consistency equations. We present numerical solu- 
tions for all temperatures including T — 0. The spin-static 
T = critical point well agrees with the results for an ear- 
lier model version 

The main results follow in Sec. We evaluate the 
critical line in the T—t plane (t represents the hopping 
strength) as a sequence of improved dynamical approxi- 
mations, which helps to derive the decay of the critical 
temperature towards the quantum critical point. This dy- 
namical approximation scheme is not designed to capture 
the zero-temperature limit. But, by increasing the number 
of Fourier components of the effective random medium, 
one finds the T c -deviation from the spin-static approxi- 
mation result almost accurately down to lower and lower 
temperatures. It is seen that the higher-frequency correc- 
tions to T c are small; their infinite number accumulates 
and leads to the non-analytical behavior of the T c -curve 
for T c — ► 0. From the characteristic decay of these correc- 
tions we deduce the location of the QCP. The numerical 
results also fit with the quantum-dynamical shift expo- 
nent of the T c -curve and thus provide a reliable estimate 
of the QCP's position. 



2 Effective action and construction of the 
self-consistency method 

2.1 Model and spin glass decoupling procedure 

We consider the grand canonical Hamiltonian 

with the fermionic Ising spin operators given by Sf = 
a ll a rf ~ a li a il- The sum index (ij) in the hopping term 
denotes summation over nearest-neighbor lattice sites. We 
assume quenched Gaussian disorder among the magnetic 
coupling constants Jy according to the distribution 

P<J «> = 7^r xp (-i ) • (2) 

In distinction to previous work [1113114) there is no dis- 
order in the kinetic part of (TjQ). In order to facilitate the 
solution of this model we assume a fully connected mag- 
netic interaction among the N lattice sites. The hopping 
takes place on an underlying simple cubic lattice in the 
limit of infinite spatial dimensions d. To obtain physically 
meaningful results we apply the usual scaling of the model 
parameters f^H] 

J = J/Vn and t = t/Vd. (3) 

To some extent the derivation of the self-consistency 
equations of the present model follows the detailed discus- 
sion given in jll| . In the following we restrict ourselves to 
a brief outline of the two-step decoupling procedure in the 
replica formalism. 

We start from the Grassmann field theoretic formula- 
tion of the n-fold replicated partition function. The dis- 
order average, i.e. integration over the magnetic coupling 
constants with the Gaussian weight J2J generates four- 
spin products in the effective action. Due to the assumed 
complete connectivity of the magnetic interaction these 
four-spin products can be reduced to quadratic terms by 
means of site-global, replica and imaginary-time depen- 
dent real decoupling fields Q r Ja ■ The disorder averaged 
replicated partition function at this stage reads 

[Z"}j = c" JvQ Jv<Pe- A °- At - A > (4) 

with the action terms 

A =J2[^La(9r-^La, (5) 
At = iJ2 I VLaVJaa, (6) 
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Here SJ a = J2 a ^ <T cr!f - ^ . is the time dependent Grass- 
mann representation of an Ising spin operator and the 
T-integrations extend from to the inverse temperature 
1/T. 

The further evaluation of relies on the elimination 
of the fields Q T Ja by means of a saddle point integration. 
The simplest but by no means trivial way to proceed would 
be the assumption of a replica-symmetric and static (i.e. 
rr'-independent ) saddle point |11I3| . Since we want to 
explore however the role played by the quantum dynamics, 
we have to keep the time dependence. 

Let us symbolize the quantum statistical average and 
the disorder average by () and respectively. Then, 
the saddle point matrix can be expressed in terms of the 
corresponding averaged spin products: 



~\TT 

I a a 



s.p. 
s.p. 



5T PT 



(Sf a )(S^)] J (8) 
■ (9) 



qt — t' cO 



Clearly, the inter-replica spin correlations are indepen- 
dent of time because the fermions can not propagate be- 
tween different replications of the system. All quantum- 
dynamical behavior of the model originates from the di- 
agonal elements of the saddle point matrix. 

In this publication we choose a global replica-symmetric 
saddle point (which is approximate only below T c , but 
does not affect the T c -result of the second order phase 
transition) with the appropriate time dependence accord- 
ing to PEJ, 

QZte =q and Q T J a = q\ T -r'\- (10) 
^ H s.p. s.p. 1 1 

We prefer to work in frequency space and perform 
Fourier transformations of the Grassmann variables as 
well as of the diagonal part of the saddle point matrix 
which take the form 



Introduction of a space saving definition. Throughout 
this article we will use the shorthand notation of the Gaus- 
sian integral operator 

£ f(x) = -j= jH dx exp (-^j f{x) (13) 

which renders many of the equations more compact. 

Employing this abbreviation we arrive at the com- 
pletely decoupled partition function 



\Z r ' 



\J 



u u n 



G 



1>1 



with the effective action 
nNJ 2 



(14) 



A 



off 



4T 2 I J2 £n -q 2 \+A t 

\ 771 / 

local 

-TV (J (H^W^ +m W^ 4- H rn *& l ~ rn \I/ 1 

/ j \ tot %<x<j iota 1 ta i a a zqo 



(15) 



In Eq. (|15fl we discarded irrelevant terms ~ n 2 . The oc- 
curring effective magnetic fields, which are complex, dy- 
namical, and local in site- and replica-indices, are given 

by 



J 



L/qzi + \/q - qyiaoj 



m>l 



'In - q.t/,.,n ) . (16) 
1 1.< ,n ' ' Wiam) ■ (17) 



2.2 The dynamical CPA approach 



r = r ^ v l e - iziT , 

1 — — OQ 



Q\r-T f \ 



oo 
m— — oo 



(ii) 

(12) 



where zi and oj m denote fermionic and bosonic Matsub- 
ara frequencies, respectively. The Fourier coefficients q m 

are real quantities and obey the symmetry relation q m = 

~ ~* 

1— m — 9 m - 

In order to facilitate the Grassmann integration the 
effective action must be further reduced to a quadratic 
form of the Grassmann variables which requires another 
decoupling step. Without going into detail we hence in- 
troduce Gaussian integrations over site-local and replica- 
global fields Zi and site-local and replica-global fields yt a . 
The latter again split into yi a o, which are also present 
in the spin-static theory, and dynamical decoupling fields 
yf am>1 , which comprise the quantum dynamical character 
of the model. 



According to Eqs. I|14l I15fl the problem has been reduced 
to an ensemble of non-interacting fermions moving in a 
complex replica- and spin-dependent effective random me- 
dium. This situation immediately calls for a dynamical 
version of the CPA [5] . Following the prescription of this 
method we replace the complex random medium iJ™ by 
a yet unknown self energy Si . The chosen limit of infinite 
spatial dimensions simplifies the problem to a single site 
problem and justifies the assumption of a site-diagonal (or 
independent) self-energy |6llj . 
The effective action (|15|l is not diagonal in the en- 
ergy indices thus allowing for virtual absorption and emis- 
sion of dynamical field quanta H™- 1 . However, the full 
fermion Green's function of the original interacting prob- 
lem with any realization of the quenched disorder is cer- 
tainly energy conserving and so is the full disorder av- 
eraged Green's function. Hence its off-diagonal elements 
in frequency space must vanish due to the average over 
the effective random medium. There is also no ferromag- 
netic tendency in the system which altogether justifies our 
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Ansatz of a spin-independent and frequency-diagonal self- 
energy. 

In order to determine Si self-consistently we keep the 
random medium at one single site, say i = 0. The effective 
action then reads 



a cpa 



y2 Z-^/ 

ra 

-tE«L(( g_1 ) 



(18) 



with the inverse of the full disorder averaged Green's func- 
tion 

(G" 1 )^ = ((izi % + tS m ) 5m (19) 

and the effective dynamical potential at the special lattice 
site i = 



V; 



i'i 



V = I + to, to > 



aH^ + Si, l' = l 



crH 7 n a *, I' = l-m,m > 0. 



(20) 



In the assumed case of infinite spatial dimensions near- 
est-neighbor hopping of non-interacting particles is de- 
scribed by the function (recall the scaling iffll) 

T (x,i) = -sign(Imx)— exp ( J (« + ern ^J 



which reveals for x = e + i0 + the well-known Gaussian 
density of states [lftllfij . Hence the full Green's function 
G defined in lfH?|) is readily expressed by 



(G)„ = T<y, 



(22) 



where T is site independent and diagonal in the Matsubara 
indices: 

(T) in =T {izi+n- E u t)6 in . (23) 

For the formulation of the self-consistency equations it 
is useful to define as an auxiliary quantity the site-local 
propagator at site i = in the presence of the effective 
potential, 



VI 

(24) 

In essence, within the CPA method the frequency-depend- 
ent self-energy Si is determined by the demand that the 
average of the local propagator matrix V aa at the special 
site i = with respect to the dynamical potential H20|l co- 
incides with the local part of the full disorder averaged ho- 
mogeneous Green's function, T. Note that equating both 
quantities requires the replica limit lim„^o to be taken. 

Omitting the algebraic details we obtain the condi- 
tional matrix equation (to condense the notation all su- 
perfluous site-, replica-, or spin-indices are dropped from 
now on) 

T = f [r4 z) (25) 

J z 




Fig. 1. Structure of the self-consistency problem of the itin- 
erant spin glass model (inspired by [Sj). (a) indicates the ef- 
fective dynamical random medium H™, which depends on the 
spin glass order parameter q and the replica-diagonal saddle 
point q m (Eqs. 1161 1171 ). (b) and (c) illustrate the determina- 
tion of the homogeneous self-energy Si by the CPA-equation 
12511 . From the local propagator f a at the special lattice site 
i = the physical quantities q and q m are constructed, which 
in turn generate the original random medium in (a). 



where the symbol [] y is a shorthand notation for the 
average with respect to all replica-local decoupling fields 

y = {VtoVmiVm}' 



w _ 0W /J, ly- W(z,y)f(z,y) 

[j \ z i y)\y 



Jy U m >lJy+Jy-W(z,y) 



(26) 



The weight function W(z, y) that appears in (|26[1 results 
from the integration of the Grassmann fields and is given 
by (Trjo- comprises spin and frequency summation) 



W(z,y) = exp(Tr i(T ln[l+V CT T]). 



(27) 



Note that, as discussed above, the off-diagonal elements 
of I" ^ as well as its spin dependence vanish exactly by 
integration. 

Expressions for the saddle point values (|TU)l can be 
obtained by construction of the spin products given by 
Eqs. H8I9|I in terms of the Grassmann fields at the special 
lattice site i = and application of Wick's theorem. After 
performing the Fourier transformation (|ll|l and taking the 
replica limit lim„_>o the set of self-consistency equations 
is thus completed by 



q = T 2 J ([Tr^rV 
r G 

J z ,,/ 



<7m = T- 



(28) 
(29) 



pi -fm, I pi' —rii.V pl+m,l' +m pi' \l } 

(T IT"' T 



y 
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In (|29|l the first term involves the product of sums over the 
m th super- and sub-diagonals and the second the trace of 
the matrix product of two factors V a shifted against each 
other about m elements along the diagonal. 

Note that in the limiting case of a vanishing hopping 
strength (i — > 0) Eqs. correctly recover the results 

of the non-itinerant model and with the additional 
choice of /i = —inT/2 J7| the equations further reduce to 
the SK-solution of the classical model 18 . 



2.3 General solution strategies 

Any attempt to solve the set of self-consistency equations 
I|25ll28ll2§|) faces the fundamental problem of the infinitely 
many quantities q m each of which effectuates correspond- 
ing Gaussian integrations via Eqs. flfil IT7I l2~7j) . In order to 
render the problem feasible we propose to keep only a few 
Fourier components q m with m — {0, . . . , M} and take the 
higher-frequency components to be zero, i.e. q m >M = 0. 
In turn this means that the Gaussian integrations associ- 
ated with the components q m >M become trivial. Below we 
will refer to this approximation scheme as the "dynamical 
approximation of order M" . 

Within this approximation scheme the quantum dy- 
namics is treated on energy scales ranging from luq = 
to u>m — 2ttTM. To estimate the quality of this approx- 
imation one has to compare the energy scales that are 
neglected to the hopping strength t as the model param- 
eter that generates the quantum dynamics. Thus we are 
led to 

t < cjm+i = 2ttT(M + 1) (30) 

as a simple criterion of validity of the M th order dynam- 
ical approximation. Hence, although it is neither a high 
temperature expansion nor an expansion in small t, the 
method works well especially for small t/T. In those re- 
gions in parameter space the approximation already at 
manageable low orders M excellently captures the effects 
of the quantum dynamics 

Another difficulty arises from the infinite extension of 
the matrices T, V (which becomes a band matrix with M 
sub- and super-diagonals (|2(J[1 in the dynamical approxi- 
mation of order M) and l~ in frequency space. Naturally, 
a numerical analysis requires the restriction to finite ma- 
trices of size 2(l c + 1) x 2(l c + 1), i.e. the matrices are con- 
structed in the limited frequency range z_; c _i to zi c . How- 
ever, there are also important contributions from higher 
frequencies that can not be neglected for accurate solu- 
tions. 

We overcome this problem by systematic asymptotic 
expansions of the self-consistency equations in terms of 
1/zi up to some feasible order O ((l/zi) K ). Here we ex- 
ploit the high-frequency asymptotics of the self-energy 

Sl ^ a * (<*r ^ + ° {h (2K+1) ) (3i) 

k=l 

where the expansion coefficients are easy to calculate 
averages of polynomials of the effective fields (I16II17JI . The 



initial values 
Q, Qm, S(, at, l c 



W(z,y), (27) 



1 

-I a*, Si, eqn. (25) 



adjust l c 



q, q m , eqns. (28, 29) 



Fig. 2. Basic iterative scheme |19| for the solution of the cou- 
pled self-consistency equations (explanation in the text). 



sums over the Matsubara frequencies which occur in Eqs. 
(|27l 1281 129fl can always be split up into a low-frequency 
main part and a high-frequency part which are separated 
by the cut-off index l c . While the matrix-structured main 
part has to be treated numerically the high-frequency con- 
tributions can be formulated in terms of asymptotic se- 
ries expansions of docile structure that permits analytical 
summation. 

The approximation of the high-frequency contributions 
by asymptotic series expansions introduces some error. 
The cut-off index l c has to be chosen such that this er- 
ror undershoots some given threshold of insignificance. In 
practical calculations we used different methods to de- 
termine l c . A simple way is to make trial variations of 
the matrix size at each iteration cycle and to adjust (in- 
crease or decrease) l c according to the corresponding vari- 
ations of all relevant intermediate quantities. A more di- 
rect method is to evaluate the contributions of the first 
neglected asymptotic order, i.e. O ((1/ zi) K+1 ^ , as a func- 
tion of l c and to apply some suitable smallness criterion. 
The latter method turned out to be un-practical for M > 
because of the complexity of the occurring analytical ex- 
pressions for the asymptotic Matsubara sums. 

The final criterion for l c is always that the physical 
quantities are independent of this auxiliary parameter at 
some desired level of precision. The proper choice of l c and 
thereupon the computational expenses for solving the self- 
consistency equations strongly depend on the temperature 
as well as on the order K up to which the asymptotic series 
expansions of the equations can be driven. 

All solutions of the self-consistency equations that are 
presented in this article have been obtained by means of 
the principal iterative algorithm sketched in Fig. [21 This 
procedure proved to be insensitive to the initial values 
and showed quite satisfying convergence properties in all 
regions of the parameter space explored so far. 

For the sake of simplicity from now on we restrict our- 
selves to the case of a vanishing chemical potential, fi — 0, 
which corresponds to half fermion filling due to particle- 
hole-symmctry of the Hamiltonian (JTJ- Without loss of 
generality we always set J = 1. 
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3 Spin-static approximation 

This section is devoted to a discussion of the dynamical 
self-consistency equations in the static approximation as 
the first and simplest of a sequence of the dynamical ap- 
proximations proposed in the last section. This static ap- 
proximation consists in neglecting the time dependence of 
the saddle point (|10[l or equivalently in taking all Fourier 
components q m with m > to be zero, i.e. M = 0. This 
restriction to the static component go implicates tremen- 
dous simplifications of the self-consistency equations. 

Because the dynamical effective fields (|17H vanish the 
decoupling fields can be integrated out trivially. One 
is left with only two Gaussian integrations over the static 
fields y and z. Also, within this approximation the occur- 
ring matrices become diagonal and thus the matrix struc- 
ture of the self-consistency equations disappears. Thus, 
the dynamical CPA-equation i|25|) decouples into a set of 
scalar equations for each Matsubara frequency that can be 
solved one at a time. Furthermore, the matrix inversion 
<|19|) turns into simple scalar inversion and the evaluation 
of the weight function (|27(l reduces to an easily manage- 
able numerical Matsubara product, 



W ata , tic (z,y Q ) = i (cosh (H /T) + 1) 



(32) 



n 

1=0 



zf + Ho 2 



where Hq — ^fqz + \Zqo~~qy - The first term in i|32|l 
is the suitable regularized frequency and spin product of 
izi + gHq and constitutes the weight function of the non- 
itinerant model In the present itinerant model the 
effect of the hopping becomes noticeable in the deviation 
of ui from izi, where the first is defined by 



ui 



T (izi - Ei,t) 
In terms of the functional 



Si 



1 



a{x) = 4x y~] — 5 — 

; = M +X 2 



(33) 



(34) 



the expressions for the saddle point values q and go given 
in (J2HI EH simplify to 



la 



T 2 £ [A(H ) 2 +A\H )] 



0) 

Va 



(35) 
(36) 



The numerical solutions of the static set of self-consist- 
ency equations are presented in Fig. [3] As T — > oo all avail- 
able many particle states become equally populated and 
go — * 1/2 since two of the four local states are magnetic. 
In the non-itinerant limit t = the paramagnet to spin 
glass transition occurs at T c = 1/ (1 + exp (— 1/(2T C ))) ~ 
0.6767. The hopping hampers the local freezing of the 
spins and lowers the critical temperature. 




Fig. 3. Finite temperature results in the spin-static approxi- 
mation for hopping strengths t = ... 1.1 in steps of 0.1. The 
upper plot shows the spin glass order parameter q (full lines) 
and the zero frequency component of the replica diagonal sad- 
dle point value go (dashed lines). The go-curves approach 1/2 
as T — > oo. Below is plotted the corresponding local static sus- 
ceptibility xo = (<?o — q)/T which remains finite as T — > 0. At 
the spin glass-paramagnet transition xo always reaches unity. 



3.1 The limit of zero temperature 

The numerical solutions of the static equations feature the 
low temperature behavior (go— g) ~ T. In order to perform 
the zero-temperature limit it is advisable to eliminate go 
and to formulate the equations in terms of g and the static 
part of the local susceptibility \o — {Qo ~ Q)/T which 
remains finite as T — > (Fig. |2J . We replace the discrete 
frequencies Zi by the continuous variable £ and recall the 
zero-temperature limit of the Matsubara summations, 



Tj2f(^i)^— dc/co- 



(37) 



By means of the rescaling of the integration variable 
Va —> Vo tne weight function (jH2J together with the 
Gaussian factor assume the form exp {—g{z, y )/T). The 
exponent function g(z,y ) can be shown to remain finite as 
T — > 0. Thus, the ^-integration reduces to a simple sad- 
dle point integration where the z-dependent saddle point 
has to be determined numerically. Finally, we derive the 
following set of self-consistency equations at zero temper- 
ature: 

To (*C- 27(0,*) = J f G E 7—— • (>!Si 

<T=±1 



u(C) + <rri(z) ' 
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0.6- 

o 
X 

<* 0.4- 



Xo 



t 

Fig. 4. Zero-temperature results for the spin glass order pa- 
rameter q — qo and the local static susceptibility xo- Increasing 
hopping strength depresses the spin glass order and drives a 
zero-temperature phase transition at t cs — 1.406 



SG, q > 




<7o 



Xo = \ — 
irq 



A'( v (z)) 



(39) 

(40) 
(41) 

(42) 



In Eqs. I|38H42|1 u(() and £((,) are continuous versions of 
the quantities ui and Si denned in Sec. 12.21 and (|33|l . In 

terms of the constant 



1-xA'Mz)) 
j](z) = xA(t](z)) + y/q\z\, 
with the functional 

A(x) = — / dC 5 . 

* Jo |«CC)I 



7 = 



^°°dCexp (2C 2 )rQ, 



(43) 



the static quantum critical point is located at t cs — 2j/tt h 
1.406. Up to a slight deviation due to different model def- 
initions this critical value quantitatively agrees with the 
results obtained in ^ for the case of a semi-elliptic free 
energy band. 

In the zero-temperature disordered phase, i.e. for t > 
t cs , the static part of the local susceptibility is given by 
Xo(t)\T=o = (Trf - \f^ 2 t 2 - 47 2 )/(27). Its deviation from 
the corresponding quantity in the non- interacting limit, 
Xo(t)\T=o,j=o — 7/( 7r 0' signalizes the vicinity of the spin 
glass phase. 



4 Quantum dynamical solutions 

While exact in the non- itinerant limit, the spin-static ap- 
proximation discussed in Sec. turns out to yield a very 
good description of our model for weak and moderate hop- 
ping. For stronger hopping, however, this static approxi- 
mation becomes increasingly inaccurate and particularly 
fails close to the T = quantum phase transition. In this 



Fig. 5. The critical line of the spin glass (SG) to paramag- 
net (PM) phase transition in the static approximation (dashed 
line) and in the dynamical approximations (full lines) of first 
(T c ,i, uppermost) to third order as discussed in Sec. 12.31 
The dotted line indicates the expected fully dynamical phase 
boundary; the light-dotted part shows the pure critical behav- 
ior l|45|l . 



section we present improved solutions of the self-consisten- 
cy equations within the dynamical approximation of up to 
third order as introduced in Sec. 12.31 

For the time being we restrict ourselves to the deter- 
mination of the phase diagram of the model. Since at crit- 
icality there is no issue of replica symmetry breaking the 
choice of the replica-symmetric saddle point (|10|) is justi- 
fied in these calculations. We determine the critical curve 
by virtue of the exact relation 



T c = qo {T c 



(44) 



which was first derived in |11| in the static approximation 
and can be shown to hold within the present dynamical 
treatment, too (see IA.lt . Equation l|44|l implies that the 
static part of the local susceptibility, xo reaches unity at 
the phase transition. In solving the conditional equation 
(|44fl it is sufficient to fix the spin glass order parameter 
to q = thus rendering the z-integrations in Eqs. I|25II29|I 
trivial. 

Our solutions for the critical line T c (t) as a sequence 
of the first three orders of the dynamical approximation 
(M = {1,2,3}) are shown in Fig. [SJ With increasing hop- 
ping strength and decreasing temperature the growing in- 
fluence of the discrete dynamic saddle point components 
Qm>o is getting more and more apparent. It can be seen 
clearly from Fig. H3 that with increasing order of the dy- 
namical approximation two successive solutions start to 
separate at larger t. We observe a rapid convergence of this 
sequence of solutions except for the region where the quan- 
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Fig. 6. Differences of the critical lines in two successive orders 
of the dynamical approximation, Am ~ T c> m-i — Tc,m- 



turn phase transition is expected. As T c — > all curves 
collapse into the static critical point at t cs . 

In the disordered phase the q m and consequently the 
effective potential matrix 1|20[1 vanish linearly with tem- 
perature, i.e. the dynamical susceptibility \ m = q rn /T has 
a finite zero-temperature limit. Hence the self-consistent 
inclusion of any finite number of the q m can affect neither 
the position of the QCP nor the critical exponents. In or- 
der to capture the quantum dynamical character of the 
problem it is necessary to take into account the q m over a 
finite range of Matsubara frequencies u> m around u>o = 0. 
Thus, the linear temperature decrease of the q m is com- 
pensated by the increasing number of Fourier components 
within a fixed frequency range and the location of the 
critical point is shifted towards smaller hopping strength 
compared to the static approximation fSec. IA.2|) . 

Close to zero temperature the critical line behaves like 

T c ~(t c -i)*, t<t c (45) 

with the shift exponent changing from <f> = 1 in the static 
approximation to tfi = 2/3 0] due to the quantum dy- 
namics (Sec. IA.2|I . Figure El gives an impression of how 
this non-analytical behavior emerges from the sequence 
of the (analytical) approximate solutions. The differences 
between two successive approximations, 

Am = T C) m-i — T C: m, (46) 

exhibit pronounced maxima (Fig.^J. While the positions 
of these maxima vary only very little they become lower 
in height but sharper with increasing M. 

The critical lines T Cj m are monotonically decreasing 
functions of t. Hence the distance between T Cj m and the 
fully dynamical true critical line, A M = T c .m — Tc.oo, pos- 
sesses a non- analytical maximum exactly at t c for any M. 
Since 

oo 

A M = J2 A **'> ( 4? ) 

M'=M+1 

this non-analyticity must coincide with the position of the 
maxima of the Am as M ^ oo. This simply means that 
the sequence of the critical lines T C: m converges slowest 



in the very proximity of the QCP. Based on this scenario 
we estimate the location of the QCP: we plot the maxima 
positions vs. their heights and extrapolate to zero height. 
This simple procedure yields the final result 

t c ~ 1.30. (48) 



5 Summary and outlook 

We considered a fermionic spin glass model including a 
nearest-neighbor hopping term. By means of standard de- 
coupling techniques the problem has been reduced to that 
of a dynamical random field system. A set of self-consisten- 
cy equations for the spin glass order parameter q and the 
Fourier components of the replica-diagonal saddle point 
qm has been derived by virtue of a dynamical CPA method. 
In order to facilitate numerical solutions we kept only a 
manageable number of low frequency components q m <M 
in the equations and abstain from the self-consistent eval- 
uation of the q m >M- We referred to this scheme as the 
dynamical approximation of order M. 

Within the static approximation (q m >o = 0) we pre- 
sented solutions both at finite and zero temperature. The 
second order T = phase transition was found at t cs = 
~ 1.406, Eq. H43fl . In order to determine the phase 
diagram of the model we calculated the SG-PM phase 
boundary in the dynamical approximation in up to third 
order (Fig.[SJ|. These data allowed to estimate the location 
of the fully dynamical critical point at t c ~ 1.30. In order 
to confirm this result it would be desirable to find solu- 
tions in higher orders of the dynamical approximation, i.e. 
M > 3. 

In this article we concentrated on the spin sector of 
our model and left out the properties in the charge sector 
such as the fermionic density of states. In the future it will 
be of high interest to investigate the effect of the hopping 
on the band structure of the system, particularly on the 
spin glass gap at zero temperature |2()j . In this context an 
extension of the solutions to non-zero chemical potential 
H is also desirable, too. 

There are also important questions concerning the in- 
terplay between quantum dynamics and replica symmetry 
breaking both issues being most significant at T = 0. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft under research project Op28/5-2 and by the 
SFB410. One of us (M. B.) also wishes to acknowledge the 
scholarship granted by the University of Wurzburg. 



Appendix 

A.l Derivation of Eq. l|33j> 

In order to locate the spin glass phase transition we ex- 
pand Eq. H28fl in terms of q (note that q enters the equa- 
tions only in the combination ^fqz by Eq. (|16|) ). With the 
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abbreviation A{z,y) = TviauV a (z, y) we have 



q = T J z \ [A] * 



q=0 



(49) 



Recall the definition of the average []y , Eq. (|26ll . The 
symmetry relations A(Q,y) — — A(0, — y) and W(0, y) = 
W(0, -y) readily follow from Eqs. (SI El EDI E3) ■ Hence, 
the first term in Eq. (|49|l vanishes by the y-integration. 
For the second term in Eq. l|49Jl we need to evaluate 



d^ z WA 



(d^ z W)A 



Wd^ z A. 



(50) 



We expand the terms In [1 + VrT] and [T^ 1 + Vr] that 
occur in the expressions for W and A, respectively (Eqs. 
(23E3J), in powers of the matrix Vr (201 • After taking the 
derivative these series can be re-summed easily yielding 



djgzW = WA, 



O^qzA 



(51) 
(52) 



Altogether, very close to the phase transition where, q ~ 0, 
Eq. reads 



T 



(Tr^oT^) - Tri a rl 



(53) 



The remaining Gaussian z-integration evaluates to 1. By 
comparing to Eq. (|29[) at T — T c , the right hand side of 
Eq. <|53|l can be identified with eft /T c 2 which finally proves 
relation (32). 



A. 2 Expansion in small q m 

In order to extract the behavior of the critical line T c [t) at 
the quantum phase transition we expand the self-consisten- 
cy equations in powers of the effective potential matrix 
V a (25) in the disordered phase, i.e. for q = 0. The inter- 
nal summations due to the occurring matrix multiplica- 
tions compensate the linear temperature decrease of the 
Fourier components of the dynamic saddle point, q m . An 
expansion up to second order yields the following self- 
consistency equation for the local dynamical susceptibility 
Xm = q"\ m \/T: 



Xm =C m (l + Xm) 

- (1 - C m Xm)T 



(54) 



OO 



^ 1 

n— — oo 



CnXr 



Using the abbreviation Xi = 1/ (l/To(izi — Si,t) + Si) 
(see (2D) the coefficients that occur in Eq. Ip^E)) are defined 

by 

OO 

C m = -2T XiX l+m , (55) 

I — — OO 

oo 

), (56) 

I — — OO 



and the self-energy is given explicitly by 



s, = -x?t V -- 

1 2 1 ^ 1 



(57) 



By numerical evaluation of Eq. (|54|l in the zero-temperature 
limit we determine a critical hopping strength t c ~ 1.372 
well between the static result and the expected fully dy- 
namical value given in Sec. 0] 

Expansion of Eqs. (|54H5fi|) in the hopping strength 
around t c = t c (T = 0) and in small frequencies Lo m finally 
yields the condition for the critical line at small tempera- 
tures 



t c - t c (T) ~T 



T E 



.=0 



1=0 



j.3/2 



(58) 
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